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ABSTRACT 


Since  mid-1961  the  National  Bureau  of 
Standards  has  used  two   rubidium  vapor  frequency- 
standards  as  transfer  or  Working  Frequency  Standards 
(WFS).      The  WFS  provides  a  continuously  available  fre- 
quency which  is  periodically  calibrated  in  terms  of  the 
United  States  Frequency  Standard  (USFS).      This  paper 
evaluates  these  calibration  data  over  a  44.  5-month 
period  of  time.     It  gives  the  long-term  performance 
of  the  rubidium  standards  in  terms  of  the  USFS,    com- 
pares them  with  a  commercial  cesium-beam  frequency 
standard,    and  presents  their  reliability  characteristics 
as  shown  by  the  mean-time-between-failures  (MTBF). 
The  calibration  data  are  graphed  with  applicable  toler- 
ance limits  for  the  between-adjustment  periods.     The 
average  standard  error  of  estimate  about  least  squares 
lines  fitted  to  the  frequency  data  is  about  2  parts  in 
10      .      The  rubidium  standards  also  show  a  rather  con- 
sistent  daily  aging  rate  of  parts  in  10x    .     This  would 
indicate  that  the  frequency  of  either  unit  would  change 
a  few  parts  in  10    u   if  operated  continuously  for  one 
year  without  any  frequency  adjustments. 


Key  Words:     aging,    calibration  (USFS),   frequency 
standards,    gas  cells,   performance 
(long-term),    reliability,    rubidium, 
stability 
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THE  LONG-TERM  PERFORMANCE  OF  TWO 
RUBIDIUM  VAPOR  FREQUENCY  STANDARDS 

B.    E.    Blair  and  A.    H.    Morgan 


1.     INTRODUCTION 

In  the  field  of  atomic  frequency  standards  the  arrival  of  commercial 
rubidium  (Rb)  gas  cells  has  aroused  considerable  interest  in  what  one  can 
expect  of  their  long-term  performance  and  stability.      However,    as  pointed 
out  by  Farmer  LI 963 J,   there  is  a  lack  of  such  data  simply  because  of  their 
relatively  recent  appearance.     Since  the  first  two  commercially -available 
Rb  frequency  standards  have  operated  nearly  continuously  over  a  44.  5- 
month  period  at  NBS,   we  can  provide  the  record  of  their  operation  for 
appraisal  and  study.      This  paper  reports  primarily  on  such  long-term 
performance  based  on  periodic  calibrations  with  the  United  States  Fre- 
quency Standard  (USFS).     Data  on  a  third  Rb  frequency  standard,   used  for 
a  period  of  3  months  in  1964,    are  included  also. 

The  Rb  gas  cell  standards  have  been  extensively  described  in  the 
literature  [Bender,    Beaty,    and  Chi,    1958;  Whitehorn,    1959;  Andres, 
Farmer,    and  Inouye,    1959;  Carpenter,    Beaty,    Bender,    Saito,    and  Stone, 
I960;  Packard  and  Swartz,    1962;  Arditi  and  Carver,    1963;  Reder,    1963]. 
These  papers  trace  their  development  from  the  earliest  prototype  models 
to  the  present-day  commercially-constructed  units. 

This  paper  briefly  describes  the  two  particular  Rb  atomic  frequency 
standards  used  at  NBS  (called  Rb-1  and  Rb-2),  and  presents  statistics  on 
their  long-term  stability  in  terms  of  the  USFS.  It  also  includes  compar- 
ison data  between  the  Rb  standards  and  a  commercial  cesium  (Cs)  stand- 
ard for  a  6-month  period  of  time,  and  briefly  discusses  relevant  charac- 
teristics of  the  Rb  standards,  such  as  short-term  stability  and  reliability 
of  operation  as  experienced  at  NBS. 


2.     BRIEF  DESCRIPTION  OF  THE  RUBIDIUM  FREQUENCY  STANDARDS 

Early  in  1961  NBS  obtained  the  first  two  commercially -available  Rb 
atomic  frequency  standards  which,    except  for  short  periods  for  repair 
and  modification,   have  operated  continuously  from  then  until  the  present 
time.     Each  standard  occupies  nearly  3  cubic  feet  in  a  cabinet  package 
weighing  about  100  pounds.      (A  later  version  of  Rb  standard  is  now  avail- 
able in  which  the  volume,   weight,    and  power  requirements  are  all  reduced 
by  a  factor  of  1/5  or  more,    and  the  specified  performance  is  comparable 
to  the  older  model.  )     Environmentally,    the  vital  elements  of  the  instrument 
are  well  protected.     Proportionally-controlled  ovens  maintain  temperature- 
sensitive  components  within  set  tolerance  limits  of  temperature,    and  a 
triple  mu-metal  shield  surrounds  elements  susceptible  to  external  mag- 
netic fields.     A  brief  description  of  the  internal  functioning  of  these  stand- 
ards follows;  detailed  operational  features  are  available  in  the  literature 
LAndres  et  al.  ,    1959;  Packard  and  Swartz,    1962;  Arditi  and  Carver,    1963 J. 

Basically,   the  Rb  frequency  standard  consists  of  a  5-MHz  crystal 
oscillator  which  is  stabilized  by  an  atomic  resonance  of  rubidium  through 
suitable  frequency  synthesizers,   multipliers,    and  automatic  frequency 

control.     Optical  pumping  and  light -transmission  monitoring  form  the 

87 
basis  for  detecting  atomic  resonance  of  Rb       between  the  ground  state 

levels  F  =  2,   M      =  0  to  F  =  1,   M      =  0.     Figure  1   shows  a  simplified 
F  F 

block  diagram  of  the  Rb  frequency  standard.     The  heart  of  the  rubidium 

o  "7  o  c 

standard  is  the  optical  package  which  contains  an  Rb       lamp,    an  Rb 

87 
filter,    an  Rb       gas  cell  within  a  microwave  cavity,    and  a  silicon  photo 

cell  detector. 

Functional  operation  of  the  Rb  frequency  standard  is  as  follows: 

87 
after  the  light  from  the  Rb       lamp  is  directed  through  the  filter  and  gas 

cells,   it  is  detected  by  the  photo  cell.      The  intensity  of  light  transmitted 

through  the  gas  cell  to  the  detector  depends  on  the  difference  between  the 


87 
microwave  excitation  frequency  and  the  Rb       transition  frequency.     Ini- 
tially,  through  a  repetitive  process  called  optical  pumping  LKastler,    1957J, 
an  excessive  number  of  Rb  atoms  jump  from  the  F  =  1  level  to  the  upper 
F  =  2  ground  level,    accompanied  by  high  transparency  and  a  consequently 
higher  transmission  of  light  through  the  gas  cell.     Exciting  the  cavity  to 
the  critical  transition  frequency  between  the  two  ground  state  energy  levels 
(approximately  6834  MHz),    results  in  a  net  gain  of  Rb  atoms  back  to  the 
F  =  1  ground  level  and  a  degree  of  opacity  to  the  light  beam  focused  on 
the  silicon  detector.      To  provide  an  error  signal  for  correcting  the  fre- 
quency of  the  radiation  driving  the  cavity,    a  low  frequency  (107  cps)  phase 
modulation  of  the  microwave  energy  is  applied  to  the  cell  cavity.      Thus, 
the  light  reaching  the  photo  cell  fluctuates  at  this  frequency.     A  phase 
detector  compares  the  photo  cell  output  signal  with  the  input  modulation 
frequency,    and  the  resultant  phase  error  signal  continuously  locks  the 
quartz  crystal  oscillator  frequency  to  that  of  Rb  transition  through  an 

electronic  servo  system.      Thus,    all  of  the  synthesized  frequency  outputs 

87 
of  the  crystal  oscillator  are  continuously  related  to  the  Rb       transition 

frequency.     Synthesis  of  the  Rb  transition  frequency  from  the  basic  5-MHz 

signal  of  the  crystal  oscillator  occurs  through  harmonic  multiplication  and 

subtraction  using  transistors  and  varactor  diodes. 

Because  buffer  gases  are  used  in  the  Rb  vapor  cell  to  facilitate  optical 

pumping  and  to  reduce  the  Doppler  width  of  the  hyperfine  absorption  lines, 

the  effective  Rb  transition  frequency  may  be  shifted  either  up  or  down, 

depending  upon  the  type  and  pressure  of  buffer  gas  used.      Consequently, 

one  cannot  use  the  modified  transition  output  frequency  as  a  primary 

standard,   but  must  reference  it  to  a  known  frequency.     On  the  other  hand, 

with  proper  selection  of  the  buffer  gas  and  its  sealed-off  pressure  at  the 

time  of  manufacture,   one  can  obtain  any  offset  frequency  that  is  required 

-9 
for  different  time  scales.     Fine  frequency  tuning,   by  about  +5x10      ,   is 


possible  also  through  adjustment  of  the  magnetic  field.     Thus,   the  limi- 
tations which  preclude  the  use  of  the  Rb  vapor  instrument  as  a  primary- 
frequency  standard  can  prove  of  value  if  precise  matching  and  adjusting 
of  output  frequencies  is  needed. 

3.     THE  USFS  REFERENCE  BASE 

The  defined  reference  for  all  standard  frequency  measurements  is  the 

unperturbed  Cs/transition  frequency  (F  =  4,   M      =  0)  "*-*■  (F  =  3,   M      =  0) 

F  F 

of  9,  192,  631,  770  Hz  at  zero  magnetic  field  [NBS,    1964].     In  the  United 
States  this  primary  reference  frequency  has  been  derived  from  NBS-con- 
structed  Cs-beam  units  which  form  the  United  States  Frequency  Standard 
(USFS)  [Mockler  et  al.  ,    I960;  Mockler,    1964],     From  I960  to  early  1965, 


the  USFS  consisted  of  two  such  units  designated  as  NBS-1  and  NBS-2  with 

onst 
-12 


an  accuracy  of  1  part  in  10      .     Presently,   the  USFS  is  the  newly-construct- 


ed NBS-3  which  is  a  longer  Cs-beam  unit  and  accurate  to  ±1.  1x10 
(1  <t  value)  LBeehler  et  al.  ,    1966J.     These  accuracy  figures  derive  from 
the  RMS  combination  of  all  estimated  uncertainties  related  to  the  USFS 
which  may  cause  deviations  from  the  true  cesium  resonance. 

3.  1.     Function  of  Working  Frequency  Standards  and  Relation  to  USFS 

Because  the  USFS  does  not  operate  continuously,   the  NBS  established 
a  Working  Frequency  Standard  (WFS),   which  consists  of  the  Rb  standards 
and  stable  quartz-crystal  oscillators,   and  provides  standard  frequencies 
without  interruption.     Periodic  measurements  relate  the  WFS  to  the  USFS. 
From  mid-1961  to  January  1964,   the  Rb  standards  were  a  direct  link  in 
the  calibration  and  frequency  control  of  the  NBS  standard  frequency  broad- 
casts LBlair  and  Morgan,    1965],     Presently,   one  of  the  Rb  standards  in  the 
WFS  contributes  to  the  NBS-A  time  scale  [Barnes  et  al.  ,    1965;  Andrews, 
1965J  which,    since  January  1964,   has  given  time-base  control  to  such  NBS 


broadcasts.     In  addition,    the  WFS  serves  as  a  continuous  frequency  ref- 
erence for  the  inter  comparison  of  atomic  frequency  standards  via  VLF 
broadcasts  [Morgan  et  al.  ,    1965J  and  provides  an  internal  calibration 
source  to  the  NBS  Boulder  Laboratories. 

3.  2.     Methods  of  NBS  Frequency  Calibrations  in  Terms  of  USFS 

At  regular  intervals  (usually  daily  except  weekends  and  holidays)  the 
Rb  standards  were  measured  in  terms  of  the  USFS.      The  NBS  Boulder 
Laboratories  have  used  two  methods  of  frequency  calibration  with  compa- 
rable results  LBeehler  et  al.  ,    1962 J.     In  the  initially  established  procedure, 
called  the  "manual"  method,   the  5-MHz  output  of  the  standard  under  cali- 
bration directly  drives  the  multiplier  chain  of  the  USFS.     For  each  meas- 
urement,  the  operator  selects  an  upper  and  lower  Intermediate  Frequency 
(IF),    equally  offset  from  the  cesium  resonance  center  frequency,    by 
judging  the  equality  of  the  atomic  beam  detector  current  for  each  condition. 
For  this  method  each  reported  value  is  the  grand  mean  of  from  20  to  50 
independent  measurements  made  over  a  period  of  1  5  to  30  minutes.     Typi- 
cally,   such  measurements  show  a  standard  deviation  less  than  1  part  in 
10      .     (It  should  be  pointed  out  that  the  manual  method  has  not  been  used 
routinely  since  early  1963.  ) 

In  the  second  or  "servo"  method,    similar  multiplication  and  klystron 
control  circuitry  is  used  as  in  the  manual  method.     In  this  case,    however, 
the  signal  from  a  stable  transfer  oscillator  becomes  the  input  to  the  USFS 
multiplier  chain  and  its  frequency  is  then  stabilized  by  the  cesium  reso- 
nance by  means  of  an  electronic  servo  system.     Frequency  comparisons 
are  made  between  the  transfer  oscillator  and  the  Rb  frequency  standard  by 
measuring  the  average  period  of  the  difference  frequency  between  them. 
In  this  case,   the  reported  daily  value  is  the  grand  mean  of  about  ten  200- 
second  averages  taken  over  a  measurement  period  of  about  1/2  hour. 


12 
The  standard  deviation  of  such  values  is  near  2  parts  in  10       for  a  typical 

measurement  period.     (Standard  deviations  of  the  means  of  such  measure- 

ments  in  terms  of  the  USFS  are  near  5  parts  in  10        [Beehler  and  Glaze, 

1965].) 

In  a  third  measurement  methods   the  Rb  signal  is  continuously  com- 
pared with  a  drift-corrected  oscillator  (DCO)  in  a  phase-lock  servo  system 
[Looney,    1 961  J.     Presently,   the  DCO  is  the  direct  controlling  link  of  the 
NBS  low  frequency  broadcasts  from  Ft.    Collins,    Colorado  [Andrews,    1965; 
Barnes  et  al.  ,    1965J.     The  frequency  of  the  DCO  is  stable  to  several  parts 
in  10    "of  the  zero  offset  frequency  and  is  periodically  compared  to  the 
USFS.     By  this  method,   used  frequently  in  1964,   the  daily  value  of  the  Rb 
frequency  was  determined  by  integrating  the  phase  analogue  records  over 
24 -hour  periods. 

One  additional  method  of  assigning  daily  frequency  values  was  used 
infrequently  but  should  be  mentioned.     This  technique  consists  of  continu- 
ously comparing  the  Rb-1  and  Rb-2  signals  in  a  phase-lock  servo  system. 
When  one  standard  is  calibrated  then,   in  terms  of  the  USFS  over  a  given 
time  period,    a  frequency  value  of  the  other  standard  is  determined  for  the 

same  time  period  from  the  Rb-1  versus  Rb-2  phase  record.     Such  compar- 

12 
isons  are  made  typically  in  parts  in  10      ,    and  we  used  this  technique  for 

determining  Rb-2  values  for  parts  of  periods  5  to  7. 

4.     PERFORMANCE  DATA 

This  section  presents  the  long-term  performance  of  the  Rb  standards, 
gives  their  relative  stability  versus  a  commercial  cesium  standard  for  a 
6-month  period,    and  briefly  discusses  their  short-term  stability  and  power 
spectrum  characteristics. 


4.  1.     Long-Term  Performance 

Aging  effects  together  with  random  day-to-day  variability  most  often 
characterize  the  long-term  stability  of  a  frequency  standard.     Both  factors 
may  be  seen  in  the  data  plots  of  the  two  Rb  standards  for  the  period  August 
1961  to  April  1965,    given  in  figures  2  and  3.     Summary  statistics  of  the 
least  squares  lines  through  these  data  are  given  in  table  1.     (These  calcu- 
lations were  obtained  with  a  least  squares  statiscical  program  run  on  a  high 
speed  electronic  computer.     An  appendix  describes  this  program.  )    Each 
plotted  point,   as  noted  previously,   is  a  frequency  measurement  in  terms 
of  the  USFS  or  related  thereto.     However,   it  was  not  always  possible  to 
obtain  a  daily  reading  and  the  dotted  points  on  the  plots  indicate  such  miss- 
ing    data.     Since  the  use  of  any  interpolated  data  would  cause  an  unwar- 
ranted smoothing  effect  in  the  results,   the  analysis  consists  of  actual 
measurements  only.     On  the  other  hand,   published  values  were  omitted  at 
times  because  of  either  faulty  operation,    service  work,    or  scarcity  of  data 
for  a  given  subperiod.     Such  values,    although  true  and  correct,    do  not 
contribute  to  the  overall  picture  of  what  one  can  reasonably  expect  of  a 
normally  operating  rubidium  frequency  standard. 

From  period  to  period  a  random  aging  (drifting)  pattern  for  each  stand- 
ard is  apparent  in  figures  2  and  3.      This  is  more  graphically  shown  in  figure 
4  where  the  average  drift  rates  of  Rb-1  and  Rb-2  are  plotted  for  periods 
between  each  adjustment  during  1961-1965.     In  general,    Rb-1   shows  a 
persistent  negative  drift  rate,   whereas  Rb-2  exhibits  both  positive  and 
negative  rates.     During  the  first  two  or  three  years  of  operation,   major 
component  replacements  were  made,    and  it  is  not  known  what  effect  this 
may  have  had  on  the  drift  rates.     Since  we  desired  to  keep  the  frequency 
of  the  WFS  rather  close  to  the  nominal  offset,   quite  a  few  frequency  ad- 
justments also  were  made  periodically  as  is  shown  in  the  data  plots.     (The 
nominal  offset  is  an  approximation,    agreed  upon  internationally,   to  the 


difference  in  the  rate  of  occurrence  of  time  ticks  on  the  universal  scale 
(UT2)  and  seconds  pulses  on  the  atomic,    or  ephemeris,    scale  (.Hudson, 
1965J.     Such  frequency  offsets  corresponded  to  -150  x  10  in  1961, 

1964,    and  1965;  and  -130  x  10"       in  1962  and  1963.)    The  standards  show 

13 
an  average  daily  drift  rate  of  parts  in  10        during  the  period  under  study. 

The  two  curves  in  figure  4  show  some  cross  correlation,  but  this  is  not 
too  surprising  since  the  standards  are  identical  and  operate  in  the  same 
environment. 

The  variability  about  the  least  squares  lines  is  indicated  by  the  stand- 
ard error  of  estimate,    S     ,    ,    (a  measure  of  the  scatter  in  the  observed 

y/x 

points  in  the  vertical  or  y  direction  about  the  least  squares  line)  [Crow 
et  al.  ,    1 960 J,    shown  in  table  1.      (The  standard  deviation  of  the  frequency 
values,    S    ,   is  usually  larger,    since  it  includes  the  effects  of  aging.  ) 
Table  2  gives  some  yearly  statistics  of  the  performance  of  the  two  stand- 
ards for  the  44.  5-month  period.     Interestingly,    from  year  to  year  there 

seems  to  be  a  decreasing  trend  in  S     ,    .      This  results,    perhaps,    either 

°  y/x 

through  an  improvement  in  measurement  methods,    or  replacement  of  var- 
ious components  in  the  standards,    or  both.      The  average  S     ,     for  these 

12  X 

standards  is  1  2  to  13  parts  in  10      .     This  figure  includes  the  high  initial 

variation  in  1961  and  1962  caused  by  component  failures  and  their  subse- 
quent replacement.     It  is  noted  that  for  several  periods  with  substantial 

12 

number  of  days,   n,    (n  >  20),    S     .     equals  or  is  less  than  6  parts  in  10 

y/x 

The  Rb-2  standard  experienced  relatively  small  disturbance  from 
March  1964  to  April  16,    1965.     During  this  period  Rb-2  operated  continu- 
ously and  NBS  made  six  discrete  frequency  adjustments,   with  no  compo- 
nent alterations  or  replacements.     Figure  5  gives  a  plot  of  these  Rb-2 
data  with  the  frequency  adjustments  removed  for  convenience.     A  least 

squares  line  fits  the  points  quite  well;  the  average  negative  drift  rate  is 

13 
about  8  parts  in  10        and  the  standard  error  of  estimate,    S     ,    ,   is 

y/x 


0.  10  x  10         .      There  appear  to  be  small  changes  in  the  drift  rate,    how- 
ever,   during  each  of  the  six  different  periods,    as  shown  in  figure  4,    and 

this  accounts  for  a  somewhat  high  average  S     ,    ,    for  the  overall  period. 

°       y/x 

Data  for  a  later -model  standard,    Rb-64,   were  also  obtained  for  the 
period  January  3  to  March  30,    1964,    and  are  shown  in  figure  3  and  table 

1.  Generally,   the  performance  statistics  for  Rb-64  agree  well  with  those 
of  Rb-1  and  Rb-2;  interestingly  enough,    however,    Rb-64  showed  very  little 
frequency  aging  during  this  period. 

4.  2.      Comparison  of  Two  Rubidium  and  One  Cesium 
(Commercial)  Atomic  Frequency  Standards 

For  an  approximate  6-month  period  (August  1962  to  February  1963) 
a  commercial  Cs  frequency  standard  was  also  used  as  a  part  of  the  WFS. 
Frequency  calibrations  of  this  standard  were  made  and  a  plot  of  these 
values,    together  with  those  of  Rb-1  and  Rb-2,   is  shown  in  figure  6.      (For 
ease  of  comparison,    several  discrete  frequency  adjustments  were  removed 
from  the  data  of  both  Rb  standards.  )    As  may  be  seen,    during  this  period, 
Rb-2  had  an  average  positive,   drift  rate;  however,   Rb-1  and  the  Cs  stand- 
ard show  average  negative  rates.      The  oven  heater  of  the  Cs  beam  tube 
was  not  functioning  correctly  for  much  of  this  period  and  what  effect  this 
had  on  its  apparent  drift  rate  is  not  known.     During  this  comparison  period, 

the  Rb-2  data  were  the  most  variable  with  an  S     ,     of  0.  2  parts  in  10      ;  but 

y/x 

the  Rb-1  and  the  Cs  standard  data  each  show  an  S     ,     of  slightly  more  than 

10  ^/X 

0.  1  part  in  10 

4.3.     Short-Term  Stability 

Although  the  major  burden  of  this  report  is  the  long-term  perform- 
ance of  the  Rb  standards  (for  periods  of  months  to  years),    some  mention 
will  be  made  of  their  short-term  stability.      The  best  short-term  stability 
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(standard  deviation)  that  has  been  reported  was  near  1  part  in  10     •  for 

averaging  times  of  24  seconds  for  each  sample  during  an  elapsed  meas- 
urement period  of  about  13  hours  [Packard  and  Swartz,  1962J.  Packard 
and  Swartz  Ll 9 6 Z J  found  that  the  measured  stability  is  inversely  propor- 
tional to  the  averaging  time,  due  to  both  the  low  signal-to-noise  ratio  of 
the  internal  crystal  oscillator  output  signal  and  the  phase  instability  of 
the  comparison  instrumentation. 

Analogue  records  of  the  continuous  phase  difference  between  two  fre- 
quency standards  (as  given  by  a  phase-lock  servo  system  [Morgan  and 
Andrews,    1961J),   provide  some  indication  of  their  relative  stability. 
Figure  7  shows  a  sample  analogue  record  which  gives  the  phase  difference 
between  Rb-1  and  Rb-2.      The  slope  of  the  trace  indicates  the  frequency 
difference  rate  between  the  two  standards,    and  the  minor  variations  in 
the  trace  suggest  good  relative  phase  stability,    as  obtained  with  a  servo 
averaging  time  of  nearly  one  minute.     The  standard  error  of  estimate 
about  a  least  squares  line  fitted  to  points  on  this  trace  at  10-minute  in- 
tervals is  about  0.  01  usecond,    corresponding  to  an  RMS  instability  in  fre- 
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quency  contributed  by  both  standards  of  less  than  2  parts  in  10      .     Assum- 
ing that  Rb-1  and  Rb-2  contribute  equally  to  this  variation,    one  can  com- 
pute an  instability  of  about  1  part  in  1  0       for  each  standard.     Because  of 
the  relatively  long  response  time  of  the  servo  system  (servo  motor  is 
gear-reduced  by  a  factor  of  6000:1)  and  a  paper  speed  of  6  inches  per 
hour,   it  is  difficult  to  determine  short-term  stabilities  for  periods  of  60 
seconds  or  less  from  such  records. 

The  output  power  spectrum  of  a  frequency  standard  may  also  give 
some  indication  of  its  short-term  stability.     NBS  has  made  power  spec- 
trum measurements  of  a  rubidium  standard  by  comparing  a  multiplied  5- 
MHz  output  to  an  ammonia  maser  reference  signal  at  24  GHz  [Barnes, 
1964;  Barnes  and  Mockler,    1960J.     These  measurements  gave  a  clean 
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center  spectrum  with  a  bandwidth  at  the  half -power  points  of  about  4  Hz 
[Packard  and  Swartz,    1962J.      (The  comparison  instrumentation,    however, 
may  have  affected  these  measurements  to  some  extent,    and  a  more  real- 
istic bandwidth  is  probably  less  than  2  Hz.  )    Also,   there  were  low  level 
sidebands  on  the  carrier  (40  dB  below  the  24-GHz  signal)  at  ±60,    ±107, 
and  ±120  Hz.      These  sidebands  probably  originated  in  the  Rb  standards. 
Note  especially  that  the  phase  modulation  frequency,   which  provides  an 
error  signal  in  the  Rb  standard,   is   107  Hz,     In  addition,   the  noise  level 
of  the  ammonia  maser  spectrum  analyzer  is  at  least  -47  dB  below  the  24- 
GHz  carrier  [daze,    1966], 

5.     RELIABILITY  OF  Rb-1  AND  Rb-2 

Reliability  of  electronic  equipment  is  often  shown  to  follow  an  expo- 
nential frequency  distribution  LLloyd  and  Lipow,    196 Z J,    such  as 
R    =     exp  -  T/t,   where 

R  =  Reliability  or  probability  of  survival 

T  =  Time  of  projected  operation,    hrs 

t    =  MTBF  (Mean  time  between  failure  -  hrs). 

A 

Using  the  calculated  MTBF,   t,   as  an  estimator  of  t,    one  can  show 

that,    for  a  period  of  one  year  (8760  hrs),   the  Rb-2  maximum  likelihood 

A 

estimator,   R,   is  a  factor  of  5  greater  than  that  for  Rb-1.     Such  a  com- 
parison number  will  vary  with  the  time,    T,    due  to  the  exponential  factor; 
nevertheless,   it  gives  some  measure  of  their  relative  service  performance. 
Table  3  gives  the  service  performance  data  of  the  two  frequency  stand- 
ards; bar  graphs  of  their  MTBF  are  shown  in  figure  8.     Rb-2  shows  over 
twice  the  MTBF  of  Rb-1  and  about  half  the  number  of  failures.     However, 
the  principal  cause  of  failure,   namely  the  tube-type  lamp  exciter,   was 
the  same  for  both  standards.     After  the  units  were  modified  with  a 
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transistorized  lamp  exciter,    considerably  improved  service  resulted. 
For  instance,    a  trouble-free  period  of  9040  hours  followed  the  modifi- 
cation of  Rb-2.     A  time  weighted  average  for  both  standards  covering  the 
entire  period  of  operation  gives  a  MTBF  of  5570  hours. 

From  the  performance  data  over  this  44.  5-month  study  certain  sta- 
tistical inferences  can  be  made  about  this  model  of  standard  in  general. 
That  is,    one  can  derive  confidence  intervals  which,   with  a  stated  confi- 
dence,   should  enclose  the  true  MTBF.     The  following  probability  expres- 
sion gives  such  confidence  intervals: 


A 

2nt 


<     t    < 


A 

2nt 


X 


2n;  a/ 2 


X    2n;  l-a/2 


=    90%, 


where 


a  =  0.10 

P  =  Probability 

n  =  Number  of  failures 

t  =  MTBF 

A 

t  =  Computed  estimate  of  t    (Total  time  between  failures/ 

Number  of  failures) 


X 


2n;  a  /  2 


=    Percentage  point  of  chi-square  distribution  with  2n 

2 
degrees  of  freedom;  the  probability  of  \  exceeding 

CdXX 

this  is  a/2. 


(NOTE:    y    „     is  defined  as  for  a  random  sample  of  t,    LLloyd  and 

2n  t 

Lipow,    1962;   Carroll,    1962J.      From  this  relationship  confidence  inter - 
vals  for  reliability,    R,    also  can  be  found.  ) 
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The  90%  confidence  interval,    then,   for  MTBF  of  this  model  Rb  stand- 
ard,   based  on  Rb-1  and  Rb-2  service  data,   is  3610  <  t   <   9930. 

The  failures  shown  in  table  3  are  classed  as  catastrophic.     That  is, 
the  standards  completely  stopped  operation  and/or  it  was  impossible  to 
phase-lock  the  controlling  oscillator  to  the  Rb  resonance.     In  all  fairness, 
we  should  state  that  at  various  times  new  components  were  installed  in 
the  standards,   and  there  are  transitory  effects  in  the  frequency  data  be- 
fore and  after  such  installation.     It  was  also  necessary  from  time  to  time 
to  align  the  synthesizer  and  second  harmonic  amplitude.     Rb-2  showed 
progressively  increased  time  between  alignments.     From  the  middle  of 
1962  to  early  1963,   we  note  periods  of  3000  to  4000  hours  between  align- 
ments.    From  then  to  early  1965,    such  periods  increased  to  6000  and  9000 
hours. 

The  continuous  operation  of  the  Rb  frequency  standards  under  labo- 
ratory conditions  and  over  a  long  period  of  time  at  NBS  has  shown  areas 
where  improvements  were  necessary  and  required.     (The  transistor  cir- 
cuitry replacement  of  vacuum  tubes  in  the  lamp  exciter  is  one  such  in- 
stance. )    As  further  improvements  are  continuing,   it  is  not  difficult  to 
foresee  failure-free,   unattended  operation  of  Rb  standards  for  periods 
greater  than  one  year. 

6.     CONCLUSIONS 

Over  a  period  of  almost  four  years  these  two  particular  rubidium  fre- 
quency standards  have  given  quite  reliable  service  at  NBS  in  performing 
the  many  functions  of  continuously -running,   transfer  and  reference  stand- 
ards.    After  initial  malfunctions  were  corrected,   largely  through  the  use 
of  transistorized  circuitry  in  the  lamp  exciter,   improved  operation  re- 
sulted.    For  example,    after  such  modification  Rb-2  gave  over  9000  hours 
of  trouble-free  service  from  April  12,    1964,   to  April  17,    1965. 
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In  terms  of  the  USFS,    lb-1  and  Rb-2  each  showed  aging  which  aver- 

13 
ages  to  parts  in  10       per  day  lor  the  periods  under  study.     Assuming  no 

malfunctions,   this  indicates  that  if  each  frequency  standard  were  left  un- 
adjusted,   at  the  end  of  one  year's  continuous  operation  its  frequency  would 
have  changed  only  by  parts  in  10      .     If  the  performance  of  Rb-64  during 
a  3-month  period  is  a  good  sample  of  the  performance  of  improved  models, 
their  aging  would  be  considerably  less.     The  day-to-day  average  stability 
(standard  error  of  estimate)  of  Rb-1  and  Rb-2  advanced  to  about  0.  1  part 
in  10        from  the  initial  value  of  0.  2  to  0.  3  parts  in  10      .     Further  im- 
provement is  desirable  in  the  general  areas  of  aging  and  variability. 

7.     ACKNOWLEDGMENTS 

We  wish  to  thank  Roger  Beehler  and  his  staff  for  the  many  calibrations 
of  the  Rb  standards  in  terms  of  the  USFS,    and  Vincent  Heaton,    Catherine 
Barclay,    Joan  Berube,    Jean  Petersen,    and  Ursula  Palmer  for  their  con- 
tribution to  the  data  analysis  work.     We  appreciate  the  help  of  Lyle  Edlin 
in  various  aspects  of  the  data  collection.     Thanks  are  extended  also  to 
Carole  Craig  for  her  careful  typing  of  this  paper. 


14 


8.     REFERENCES 

Andres,    J.   M.  ,    D.    J.    Farmer,    and  G.    T.    Inouye  (October  1959),   Design 
studies  for  a  rubidium  gas  cell  frequency  standard,   IRE  Trans,    on 
MIL  ELECT.,   MIL-3,   pp.    178-183. 

Andrews,  D.  H.  (December  1965),  LF-VLF  frequency  and  time  services 
of  the  National  Bureau  of  Standards,  Trans.  Instr.  and  Meas.  IEEE, 
IM-14,   No.   4,   pp.    233-Z37. 

Arditi,   M.  ,   and  T.    R.    Carver  (January  1963),    The  principles  of  the 
double  resonance  method  applied  to  gas  cell  frequency  standards, 
Proc.   IEEE,    51,   pp.    190-202. 

Barnes,    J.   A.    (June  1964),   Personal  communication. 

Barnes,  J.  A.  ,  D.  H.  Andrews,  and  D.  W.  Allan  (December  1965),  The 
NBS-A  time  scale  -  its  generation  and  dissemination,  Trans.  Instr. 
and  Meas.   IEEE,   IM-14,   No.    4,   pp.    228-232. 

Barnes,    J.   A.,    and  R.    C.   Mockler  (September  I960),    The  power  spec- 
trum and  its  importance  in  precise  frequency  measurement,   IRE 
Trans,    on  Instr.  ,   1-9,   pp.    149-155. 

Beehler,   R.    E.  ,    W.    R.   Atkinson,    L.    E.    Heim,    and  C.    S.    Snider  (De- 
cember 1962),   A  comparison  of  direct  and  servo  methods  for  utilizing 
cesium  beam  resonators  as  frequency  standards,   IRE  Trans,    on  Instr., 
1-11,   No.    3  and  4,   pp.    231-238. 

Beehler,   R.    E.  ,    and  D.    J.    Glaze  (December  1965),   Private  communica- 
tion. 


15 


Beehler,   R.    E.  ,    D.    Halford,   R.    Harrach,   D.   Allan,   D.    Glaze,    C.    Snider, 
J.    Barnes,    R.    Vessot,    H.    Peters,    J.    Vanier,   L.    Cutler,    and  L. 
Bodily  (February  1966),   An  inter  comparison  of  atomic  standards, 
Proceedings  Letters,   Proc.   IEEE,    54,   No.    2,   pp.    301-302. 

Blair,    B.    E.  ,   and  A.    H.    Morgan  (July  1965),    Control  of  WWV  and  WWVH 
standard  frequency  broadcasts  by  VLF  and  LF  signals,    Radio  Sci.    J. 
Res.   NBS  69D,   No.    7,   pp.    915-928. 

Bender,   P.   L. ,    E.    C.    Beaty,   and  A.   R.    Chi  (November  1958),   Optical 

87 
detection  of  narrow  Rb       hyperfine  absorption  lines,   Phys.    Rev. 

L.    1,   No.    9,   pp.    311-313. 

Carpenter,    R.    J.  ,    E.    C.    Beaty,   P.   L.    Bender,    S.    Saito,    and  R.    O.    Stone 
(September  I960),   A  prototype  rubidium  vapor  frequency  standard, 
IRE  Trans,    on  Instr.  ,   1-9,   pp.    132-135. 

Carroll,    J.   M.    (November  1962),    Reliability:   1962,    Electronics,    35,   No. 
48,   pp.    53-76. 

Crow,    E.    L.  ,    F.   A.    Davis,    and  M.    W.    Maxfield  (I960),    Statistics  manual 
(Dover  Publications,   Inc.,   New  York,   N.    Y.  ),   pp.    147-165. 

Farmer,   D.    J.    (May  1963),    Performance  and  application  of  gas  cell  fre- 
quency standards,    Proc.    17th  Ann.    Freq.    Con.    Sym.      (U.    S.   Army 
Res.    and  Dev.   Lab.,    Ft.    Monmouth,   New  Jersey),   pp.   449-461. 

Glaze,   D.    (April  1966),   Personal  communication. 

Hudson,    G.    E.    (August  1965),   Of  time  and  the  atom,   Physics  Today, 
pp.    34-38. 

Kastler,   A.    (June  1957),   Optical  methods  of  atomic  orientation  and  of 
magnetic  resonance,    J.    Opt.    Soc.   Am.  ,   47,   pp.    460-465. 


16 


Lloyd,    D.   K.,   M.   Lipow  (1962),   Reliability:  management,   methods  and 
mathematics      (Prentice-Hall,    Englewood  Cliffs,   New  Jersey),   pp. 

190-195. 
Looney,    C.    H.  ,    Jr.    (February  1961),   A  very  low  frequency  synchronizing 

system,   Proc.   IRE,   49,   No.    2,   pp.    448-452. 
Mockler,    R.    C    (May  1964),   Atomic  frequency  and  time  standards,   Radio 

Sci.    J.    Res.   NBS  68D,   No.    5,   pp.    523-527. 
Mockler,   R.    C  ,   R.    E.    Beehler,   and  C.    S.   Snider  (September  I960), 

Atomic  beam  frequency  standards,   IRE  Trans,    on  Instr.  ,   1-9,   pp. 

120-132. 
Morgan,   A.    H.  ,    and  D.   H.   Andrews  (April  1961),   MSthodes  et  techniques 

de  controle  des  ondes  kilometriques  et  myriametriques  aux  Boulder 

Laboratories,    Comit£  Consultatif  pour  la  Definition  de  la  Seconde 

aupr^s  du  Comit^  International  des  Poids  et  Mesures,    26  session, 

Annexe  6,    (Gauthier-Villars  &  C1S,   Paris),   pp.    68-72. 
Morgan,   A.   H.  ,    E.    L.    Crow,   and  B.    E.    Blair  (July  1965),   International 

comparison  of  atomic  frequency  standards  via  VLF  radio  signals, 

Radio  Sci.   J.   Res.   NBS  69D,   No.    7,   pp.    905-914. 
Natrella,   M.    G.    (August  1963),    Experimental  statistics,   NBS  Handbook 
91,    Chapter  22  (Notes  on  statistical  computations),   pp.    22-1  to  22-4. 

NBS  (December  1964),   World  sets  atomic  definition  of  time,   NBS  Tech. 

News  Bull.,   48,   pp.    209-210. 
Ostle,    B.    (1963),   Statistics  in  research,    (Iowa  State  University  Press, 

Ames,   Iowa),    2nd  Edition,   pp.    159-206. 


17 


Packard,   M.    E.  ,   and  B.    E.    Swartz  (December  1962),    The  optically  pumped 
rubidium  vapor  frequency  standard,   IRE  Trans,    on  Instr.  ,   I- 11,   pp. 
215-223. 

Reder,  F.  H.  (May  1963),  Achievements  and  problem  areas  of  atomic  fre- 
quency control,  Proc.  17th  Ann.  Freq.  Con.  Sym.  (U.  S.  Army  Res. 
and  Dev.    Lab.  ,    Ft.   Monmouth,   New  Jersey),   pp.    329-371. 

Whitehorn,    R.   M.    (May  1959),    Gas  cell  frequency  standards  using  buffer 
gases  and  buffer  walls,   Proc.    13th  Ann.    Freq.    Con.    Sym.     (U.   S. 
Army  Res.   and  Dev.   Lab.  ,    Ft.   Monmouth,   New  Jersey),   pp.    648-654. 


18 


9.     APPENDIX 

9.  1.     Description  of  Least  Squares  Statistical  Program 
for  Use  With  High  Speed  Electronic  Computer 

The  least  squares  program  was  written  to  obtain  simple  linear  re- 
gression coefficients  of  standard  frequency  data  which  show  an  aging  or 
drift  effect  from  day  to  day.      This  program  analysis  assumes  an  independ- 
ent variable,   y,    corresponding  to  daily  frequency  values.      The  program 
requires  a  set  of  data  such  as 

Xl    yi 
X2    y2 


x.,   y.,   where  j  =   1  to  n 

and  assumes  a  straight  line  as  the  best  fit  to  these  paired  data.      By  the 
method  of  least  squares,    one  can  determine  the  sample  regression  curve 
of  y  on  x  by  minimizing  the  sum  of  squares  of  the  vertical  (or  y)  deviations 
of  each  frequency  value  from  the  average  curve  LCrow  et  al.  ,    I960  J.      The 
equation  of  this  sample  regression  line  is  y  =  a  +  bx,   where  a  is  the  inter- 
cept and  b  is  the  slope  of  the  regression  line. 

The  program  prints  out  both  the  input  data  (given  in  parts  of  1  0        to 
2  decimal  places)  and  standard  regression  coefficients  in  summary  form. 
These  regression  coefficients  are  computed  by  formulas  given  in  most  sta- 
tistical texts  [see  Crow  et  al.  ,    I960;  Ostle,    1 963 J.     The  program  accepts 
a  maximum  of  up  to  370  days  of  y  frequency  values.     Missing  y  values  are 
neglected;  however,   the  x  values  are  advanced  the  corresponding  number 
of  missing  days;  thus,    each  frequency  value  matches  the  day  or  x  variable 
it  is  reported  for.     Although  the  input  data  are  given  to  2  decimal  places, 
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double  precision  techniques  are  employed  in  the  program  whereby  18 
significant  figures  are  carried  in  the  calculations.      The  printout  gives 
half  this  number  of  figures,    and  the  final  results  are  rounded  to  be  con- 
sistent with  the  original  data.      This  follows  standard  statistical  practice 
to  insure  that  errors  of  computation  are  small  in  relation  to  the  statistical 
quantity  being  determined  LNatrella,    1 963 J. 
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FIG.  1        BASIC     SYSTEM    OF     Rb    ATOMIC     FREQUENCY     STANDARD 
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FIG.  2       CALIBRATION     DATA    OF    Rb-1     AND     Rb-2 
FROM     AUGUST    1961    TO    DECEMBER     1963 
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FIG.  3       CALIBRATION     DATA     OF   Rb-1    AND    Rb-2 
DURING    1964   TO   EARLY    1965 

(Rb-64  DATA    GIVEN   IN    PERIOD    13) 
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FIG.  6      FREQUENCY     COMPARISONS    OF    Rb    AND    Cs    ATOMIC    FREQUENCY 
STANDARDS      FOR    6  MONTH     PERIOD.     (TERMS    OF    USFS) 
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TABLE    1 
SUMMARY     STATISTICS    OF    Rb     STANDARDS     FREQUENCY     DATA 


PERIOD 

Rb  - 

1     STATISTICS 

Rb  -  2      STATISTICS 
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b 
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+  0.021 
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<  0.001 

0.22 
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* 
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rO 
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Ol 

10  a 

20 
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0  09 
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10b 
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1  1 
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0.10 
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13a 

* 

m 

r^#/ 

yA^// 
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0.96 

W/, 

/yiy 

Wd 

13b 
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0.19 
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m 

'4¥// 

w// 
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^M 

M 

W// 
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0.10 

008 
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13d 
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0.005 

0.09 
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-0.46 
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o> 

14 
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0  003 

0.07 

005 

DC. 37 

0.74 
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13 

-0.018 

0.002 

0  07 

0.02 
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15 
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STATIST 
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** 

ICS    NOTATION      (DATA    GIVEN   IN    PTS.   IN    10      —  1961,   1964    AND    1965    DATA   ADJUSTED 
TO    NOMINAL     OFFSET      OF     -130.00    X    IO~10  ) 

NUMBER    DAILY   REPORTED  CALIBRATIONS 

SLOPE    OF  LEAST   SQUARES    LINE   THRU   DATA    (DRIFT   PER    DAY) 

STANDARD    DEVIATION    OF    SLOPE 

STANDARD    DEVIATION    OF  DAILY    VALUES      (  INCLUDING  SLOPE  CONTRI BUTION  ) 

STANDARD    ERROR   OF   ESTIMATE    (  VARIATION   ABOUT   LEAST    SQUARES   LINE) 

AVERAGE    DAILY    VALUE 
^COEFFICIENT    OF    CORRELATION 
■  INTERCEPT    OF   LEAST    SQUARES  LINE 
OUT   OF  SERVICE    FOR    MODIFICATION 
Rb  -64     DATA   -  (Rb-2    IN    REPAIR  ) 


29 


CM 


x> 

rf 

13 

fl 

rt 

r— 1 

X) 

CM 

tf 

w 

<4H 

J 

o 

« 

en 

O 

■  H 

4-1 

0) 

•  H 

+J 

nj 

4-> 

CO 

>^ 

^H 

^ 

nj 

CD 

PH 

X 

1 — 1 

vO 

o 

r- 

in 

CM 

"^ 

CM 

1 — 1 

1 — 1 

o 

o 

■ — I 

>> 

a 

o 

• 

0 

, 

, 

CO 

o 

o 

o 

o 

o 

O 

00 

r>- 

OO 

v£> 

o 

00 

>^ 

CM 

CM 

. — i 

CM 

o 

CM 

CO 

o 

. 

0 

0 

, 

o 

o 

o 

o 

o 

o 

o 

■* 

1 — 1 

1 — 1 

1 — 1 

CM 

CM 

o 

o 

o 

o 

O 

o 

1 

X 

o 

o 

o 

o 

O 

o 

rQ 

CO 

o 

o 

o 

o 

o 

o 

Ph 

V 

CM 

CM 

CM 

o 

00 

00 

o 

O 

o 

1 — 1 

o 

o 

o 

O 

o 

o 

o 

o 

£1 

• 

o 

o 

O 

o 

o 

o 

o 

+ 

1 

+ 

I 

1 

1 

o 

CM 

<* 

00 

00 

CM 

S 

o 

"^ 

o 

Lfl 

^o 

lO 

CM 

i— H 

CM 

00 

X 

■tf 

i—i 

o 

00 

1 

OO 

"■ — 

CM 

■ — I 

o 

o 

i — i 

>. 

a 

* 

. 

„ 

1 

. 

CO 

o 

O 

o 

o 

o 

CM 

o 

o 

OO 

lO 

>^ 

<tf 

CM 

CM 

CM 

CM 

CO 

, 

a 

„ 

,, 

1 

„ 

O 

o 

o 

o 

O 

oO 

7—1 

CM 

<tf 

CM 

O 

o 

O 

o 

O 

> 1 

| 

X 

O 

o 

O 

o 

1 

O 

X> 

CO 

O 

o 

O 

o 

o 

ti 

V 

ro 

h- 

CM 

o 

o 

o 

o 

i — I 

1 — 1 

r— 1 

o 

o 

o 

o 

. 

o 

XI 

o 
1 

o 
1 

o 

o 

o 
1 

oo 

o 

CM 

CO 

oo 

a 

r- 

CM 

O 

o 

1 — 1 

1 

1 — 1 

J 

c<4 
< 

i — i 

CM 

0O 

■* 

LTl 

< 

xO 

vD 

sO 

sO 

O 

w 

o^ 

O 

o> 

O 

O 

>H 

1 — 1 

■ — i 

. — 1 

i — i 

■ — i 

> 

o 

<u 

h 

<D 

XJ 

= 

£ 

r— 1 

<U 

CM 

1 — 1 

CO 

nj 

H 

CO 

<u 

fl 

u 

a 

CO 

nj 

rtl 

•  H 

h 

w 

nj 

O 

> 

XI 

T3 

a 

0> 

>> 

o 

co 

o 

ft 

<+H 

) 

o 

o 

CO 

r— 1 

s 

a 

Pd 

•  H 

cu 

• 

*H 

W 

rri 

4-> 

a, 

CO 

n 

Pi 

•  H 

u 

•r-( 

n 

4-1 

ci> 

nj 

> 

•rH 
> 

OB 

T3 

o 


Ph 


c; 

T3 

•i-* 

l-i 

-u 

nl 

en 

Tj 

4-> 

s 

nj 

nj 

4-1 

4-> 

CO 

CO 

CM 


CM 


CO 


rf* 


tf* 


CM 


CM 


CO 


+ 


CM 


CM 


Pi 

+ 


CO 


CM 


Ph 


CO 


30 


co 

H 


cm 


i 

O 

CD 

u 

£ 

a 

o 

"4-1 
u 
<u 

Ph 

CD 

o 
•1-1 

> 

<d 
w 


<  H  Q 

H  <  O 

O  tf  3 

^  S?  H 

o  ft 


mi 
a  o 

H  H 

a  ^ 
o  o 


O 


w 

H 

z 

a 

n 

< 


o 


en 

H 

D 
i— i 


Q 


CO 

to 

H 

u 

43 

45 

o 

<M 

CO 

IS) 

cm 

LO 

M 

« 

lO 

CM 

N 

co 

U 

45 

o 

CO 

co 


xD 


co 
H 
45 

o 
o> 

o 

CO 


42 


N 

■ 

O 
O 


CO 

45 

o 
o 
o 

CO 


■tf 


CO 

u 
45 

o 

o 

CO 


CM 

I 
43 


co 

^ 

■y 

H 

S 

0 

<u 

+-> 

a 

0 

rd 

(5 

m 

00 

•H 

1— 1 

m 

0 

03 

cu 

CD 

a 

CU 

•rH 

4-> 

£ 

CU 

45 

43 

co 

cu 

<U 

s 

•H 

T3 

0 

• 

0) 

-t-> 

1— 1 

a 

•rH 

+-> 

CD 

0 

•r-t 

rd 

^i 

CU 

co 

CU 

> 

I— 1 

5 

rd 

+-> 

rd 

> 

co 

0) 

4^ 

u 
cu 

■H 

a 

co 

Ph 

o 

'CU 

d 

cu 

rd 

cu 

3 

U 

45 

i — i 

rd 

£ 

H 

> 

a 

u 

■w 

• 

CU 

• 

co 

cu 

(0 

Pi 

cu 
> 

0 

i-H 

CD 

a 

•H 

cu 
45 

cu 

CU 

45 

CU 

cu 

u 
rd 

rd 

■H 

rd 

co 

0 

o 

(5 

•rH 

CU 

r— I 

rd 

•H 

4-> 

nj 
U 

•  i-i 

CU 

> 

Em 

T3 

3 

T— 1 

rd 

0 

3 

* 

0 

CJ 

> 

H 

a 

45 

r-H  CM 


CO| 


31 


GPO    858-  646 


U.S.  DEPARTMENT  OF  COMMERCE 
WASHINGTON,  D.C.     20230 


POSTAGE  AND  FEES  PAID 
U.S.   DEPARTMENT  OP  COMMEF 


OFFICIAL  BUSINESS 


